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Abstract We study the effects of outflow/wind on the 
gravitational stability of accretion discs around super- 
massive black holes using a set of analytical steady- 
state solutions. Mass-loss rate by the outflow from the 
disc is assumed to be a power-law of the radial distance 
and the amount of the energy and the angular momen- 
tum which are carried away by the wind are parame- 
terized phenomenologically. We show that the mass of 
the first clumps at the self-gravitating radius linearly 
decreases with the total mass-loss rate of the outflow. 
Except for the case of small viscosity and high accretion 
rate, generally, the self-gravitating radius increases as 
the amount of mass-loss by the outflow increases. Our 
solutions show that as more angular momentum is lost 
by the outflow, then reduction to the mass of the first 
clumps is more significant. 

Keywords galaxies: active - black hole: physics - ac- 
cretion discs 

1 Introduction 

Accretion discs are believed to be present in a wide va- 
riety of astronomical systems. The winds/outflows are 
potentially important as they may enhance the accre- 
tion rate through carrying away angular momentum. 
It is now widely accepted that winds or outflows have 
their origin in accreting systems, which at same time 
power the radiation emission associated with the ob- 
ject (e.g., Blandford & Payne 1982; Fender, Bclloni & 
Gallo 2004). In a few cases, outflows have been demon- 
stratively observed in emission both from the broad line 
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regime in narrow-line Seyfert 1 galaxies (e.g., Lcighly 
& Moore 2004; Leighly 2004) and from the narrow line 
region of Seyfert 1 galaxies (e.g., Das et al 2005; Das et 
al 2006). Outflows from Active Galactic Nuclei (AGN) 
are much more common than previously thought: the 
overal fraction of AGNs with outflows is fairly constant, 
approximately 60%, over many order of magnitude in 
luminosity (Ganguly & Brotherton 2008). 

The possible effects of winds/outflows on the ra- 
dial structure of an accretion disc can be studied 
(e.g., Knigge 1999; Combet & Ferreira 2008) within 
the framework of standard theory of accretion discs 
(Shakura & Sunyaev 1973). It seems that outflows are 
more significant in the outer parts of a disc. But accord- 
ing to the standard theory of accretion discs (Shakura 
& Sunyaev 1973), the outer parts of steady, geometri- 
cally thin, optically thick discs are prone to self-gravity 
and they might be expected to fragment into stars (e.g., 
Shlosman & Begelman 1987; Goodman 2003). While 
many authors proposed possible solutions to this prob- 
lem (e.g., Goodman 2003), it seems such a gravitation- 
ally unstable disc is a good explanation for the existence 
of the first starts in galactic centers or the formation of 
super massive stars in quasar accretion discs (Good- 
man & Tan 2004; Tan & Blackman 2005; Levin 2007). 
However, as far as we know, the effects of outflow on 
the gravitational stability of accretion discs of ANGs or 
around supermassive black holes have been neglected 
for simplicity. Many authors have studied the gravita- 
tional stability of discs without outflows during recent 
years. For example, Goodman & Tan (2003) estimated 
the mass of first clumps in quasar discs. Assuming that 
the viscosity is proportional to gas pressure, they found 
that the self-gravitating radius is ranging from 1700 to 
2700 Schwarzschild radius, with a few hundred solar 
mass for the fragments. In another study, Nayakshin 
(2006) studied star formation near to our Galactic cen- 
ter. The mass of the first stars of this model is a few so- 
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lar mass which may increase because of the subsequent 
accretion. Also, Khajenabi & Shadmehri (2007) studied 
accretion discs around quasi-stellar objects (QSOs) and 
the Galactic center with a corona and determined the 
self-gravitating radius and the mass of the first clumps. 
They showed that existence of a corona implies a more 
gravitationally unstable disc. 

The purpose of the present work is to study possi- 
ble effects of outflows on the gravitational stability of 
accretion discs around supermassive black holes. We 
present a set of steady-state analytical solution for the 
structure of an accretion disc with outflow. Basic as- 
sumptions and the main equations are presented in the 
next section. Gravitational stability of the disc is an- 
alyzed in section 3 by determining the self-gravitating 
radius and the mass of the first clumps. Summary of 
the astrophysical implications is given in section 4. 



2 General Formulation 

We are following an approach similar to Knigge (1999) 
in order to include outflows in the main equations. We 
can start by defining the cumulative mass-loss rate from 
the disc as 
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where i?; n denotes radius at the inner edge of the disc 
and m w (i?) is the mass- lose rate per unit area from 
each disc face. Having this definition, the total mass 
loss rate is given by 



M w ,totai(i?) - M w (R diac ), (2) 

where i?disc is the radius of the accretion disc. 
The equation of mass conservation is written as 
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where Vr is the radial inflow velocity of material in the 
disc (Vr < 0). Note that the last term on right hand 
side of equation of mass conservation represents mass 
loss by the outflow. 

Also, we can write equation of conservation of angu- 
lar momentum. 
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The first two terms on the right-hand side of this 
equation describe the inflow of the angular momentum 
through the boundaries of the annulus and the effects 
of viscous torques due to shear (here, T r ^ is viscous 
stress). The third term allows for the angular momen- 
tum sink provided by the outflow. Parameter I permits 
the most meaningful division of the available parameter 
space and, in principle, allows different types of accre- 
tion disc winds to be identified. Choice of I = 1 should 
be most appropriate for radiation-driven outflows (e.g. 
Murray & Chiang 1996; Proga et al. 1998) and I > 1 
is appropriate for ccntrifugally driven MHD winds (e.g. 
Blandford & Payne 1982; Pelletier & Pudritz 1992). 
Vertical hydrostatic equilibrium of the disc implies 
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and to keep the model simple, the disc will now be 
assumed to be Keplerian and in a steady state. The 
first of these assumptions gives 



GM , 
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where M is the mass of the central object, whereas 
the second implies that d/dt = in ([3]) and J3J. The 
conservation equation §3§ can then be integrated as 



-2tvKEVr - M W (R) = constant. 



(7) 



The first term on the left-hand side of this expression is 
clearly the accretion rate at radius R, i.e. M acc (i?). The 
constant of integration in equation J7]) can be found by 
considering the boundary condition at R = i?j n , which 
shows that it is simply the rate of accretion on to the 
central object, i.e. 



M acc (i?) = M acc (i? in ) + M W {R). 



(8) 



The angular momentum integral is somewhat more dif- 
ficult to compute, since a result of the factor R 2 fl in 
the new outflow sink term, the right-hand side of equa- 
tion ((D) is no longer a perfect differential (Knigge 1999). 
That term , at least, must therefore be integrated ex- 
plicitly, which requires M W (R), or equivalently, m w (i?) 
to be specified. 

A simple power-law profile has been used by many 
authors (e.g., Quataert & Narayan 1999; Knigge 1999; 
Beckert 2000; Turolla & Dullemond 2000; Misra & 
Taam 2001; Fukue 2004). So, 



,{R) = KRt, 
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where the constant K being fixed by normalizing to 
M w , total (which is used as a free parameter) and £ that 
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is the radial mass-loss power-law index which is used as 
a free parameter also. 

A more general prescription for the viscous stress T r 
is considering (Taam & Lin 1984; Watarai & Mineshige 
2003; Merloni & Nayakshin 2006): 



l-u/2 a/2 
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where ao and < fi < 2 are constants and p is the sum 
of the gas and radiation pressures. Also, /3 is the ratio 
of gas pressure to the total pressure. 

Having equation (|10| as prescription for the viscous 
stress and using simple power- law form for m w (R), 
equation ([5]), then equation (fj| can be integrated 
straightforwardly. After mathematical manipulations, 
we can write equation ((¥]) as 
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The energy balance of the disc is given by 
D(R)=aT^ + e w {R), 



(12) 



where D{R) is the viscous dissipation in the disc and 
e w (i?) represents energy- loss rate by the outflow. As in 



Knigge (1999), we parameterize the outflow loss term 
e w : the rate at which energy must be supplied to out- 
flow in order for it to overcome its binding energy, i.e 
e w = (l/2)ijbrii w i;^, where r/^ < 1 is a efficiency con- 
stant. Also, kinetic energy has to be supplied to the 
outflow at a rate of ek = (l/2)r]] i rh w v' 2 a , where % < 1 
is a efficiency constant and Voo(R) = /vk (/ is a con- 
stant). Thus, equation (fT2"|) becomes 

D(R) = aT c 4 ff + i(r ?b + ^/Vw^OR), (13) 

and finally after mathematical manipulation, the en- 
ergy equation (|13[1 is written as 
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Finally, with the vertical transport of heat domi- 
nated by radiative diffusion we have a relation between 
the midplane and surface temperatures which is given 

by 



T = (^S) 1 /4 Tcff , 



(16) 



where k is the opacity coefficient. 

Equations ©, (QT|), (O and (JTHJ) enable us to find 
p and T and p as functions of R with (3 the critical 
input parameters. Thus, 
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(XpckB \— 2/3 and the ratio j3 is obtained from nondimensional form 

P ~ 3na 3cr^ m m H A4 cc (i? disc ) 2 r24 j2 f equation ((20]), i.e. 

x (1 - /3 )-2/3 /3 (8-5^)/i 2(1 _ ^ )5 /6 (1 _ $ _ e) i/ 2j (lg) 0.16^ /6 /i 3 / 2 Af 8 1/6 (^) 4 / 3 J 4 / 3 r 3 " 7/4 



87rao 



" 3f2^ Af acc (-Rdisc )J(R) 



(19) 



There is an algebraic equation for /3 as follows 
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where ^ m is the mean particle mass in units of the 
hydrogen atom mass, mu- The other constants have 
their usual meanings. In order to study the behavior 
of our solutions, it is more convenient to introduce di- 
mensionless variables. For the central mass M, we in- 
troduce Mg = M/(10 8 Mq) and for the radial distance 
R, we have r 3 = R/(1Q 3 R S ), where R s = 2GM/c 2 is 
the Schwarzschild radius. The mass accretion rate can 
be written as 
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where Ze = L/Le is the dimensionless disc luminosity 
relative to the Eddington limit, e = L/(Mc 2 ) is the 
radiative efficiency and K e . s . ~ 0-04 m 2 kg~ 1 is the elec- 
tron opacity. 

Our solution for T, P and p are written as 
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where k = K/n e . s , We can also calculate the surface 
density as 

£ = 1.27 x 10 V /3 £Af 8 1/3 (^) 5 /3 J 5 / 3 r 3 - 11/4 

I Xj g(4-f0/3(i _ /3)- 4 /3(l - *) 2 /3(! _ $ _ e) (27) 

Toomre (1964) showed that a rotating disc is subject 
to gravitational instabilities when the Q-paramctcr be- 
comes smaller than a critical value, which is close to 
unity, 
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where c s is the sound speed inside the accretion disc 
and Q = J7k is the angular velocity. So, the Toomre 
parameter of our model becomes 

Q = 44^ /2 K- 3 / 2 Af- 3/2 (^)- 2 J- 2 r 3/2 



x/3 -(8-M)/4( 1 _ /? )2( 1 _ ^)-l/2( 1 _ ^ _ 0)"3/2_ ( 2Q ) 



This equation with algebraic equation ([26]) gives the 
Toomre parameter as a function of the radial distance. 
Note that in the case of no-outflow our solutions reduce 
to standard disc solutions. Generally, the Toomre pa- 
rameter is much higher than unity in the inner parts of 
the disc which implies these regions are gravitationally 
stable and do not fragment (e.g., Khajenabi & Shad- 
mehri 2007). But the Toomre parameter decreases with 
increasing radial distance so that Q reaches the criti- 
cal value of unity at a self-gravitating radius which we 
denote by R sg . Thus, all regions with R > R sg are 
gravitational unstable and may fragment to clumps and 
cores. 

Different authors estimate the mass of fragments dif- 
ferently. Since the disc is marginally unstable, the ini- 
tial sizes and masses of gravitationally bound fragments 
can be determined by Toomre's dynamical instability 
(Toomre 1964). The most unstable wavelength for the 
Q ~ 1 disc is of order of the disc vertical scale height 
H (Toomre 1964). Thus, the most unstable mode has 
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Fig. 1 Mass of the first clumps, Mf rag , (in solar mass) at the self-gravitating radius versus the ratio of the total mass- loss 
by the outflow to the accretion rate, M Wi totai/M acc (i?disc), for fi — 2.0, /_t m = 0.6, k = 1, Ms = 1, / = \/2, »7b = 1 and 
77k = 1- Solid and dashed curves are corresponding to the solutions with I — 1 and I = 5, respectively. Each curve is marked 
by a pair The other input parameter are a = 0.03 and fe/e = 1 (top, left), a = 0.03 and fe/e = 10 (top, right), 

q = 0.3 and Zg/e = 1 (bottom, left), a = 0.3 and Ze/e = 10 (bottom, right). The general treatment of the solutions shows 
that as more mass is extracted from the disc by outflow, i.e. a stronger outflow, the first clumps are forming with smaller 
mass. 
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Fig. 2 The same as Figure [T] but it shows self-gravitating radius, i? sg , (in Schwarzschild radius Rs) versus the ratio of the 
total mass-loss by the outflow to the accretion rate, M w ,totai/-Macc(,Rdisc)- Except for the case of small viscosity coefficient 
and high accretion rate, i.e. a = 0.03 and Ze/e = 10, generally, self-gravitating radius increases as the outflow becomes 
stronger. 
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radial wave number A: mu = (QH) 1 and so the mass of 
a fragment at R = R sg becomes 

M frag a S(-t^) 2 = 47r 2 Si7 2 . (30) 

In the next section, we calculate the self-gravitating 
radius and the mass of the first clumps according to 
our analytical solutions. 

3 Analysis 

We fix the central mass, the opacity and the mean 
molecular weight, respectively as Mg, = 1, k = 1 and 
Mm = 0.6. Also we have / = y/2, 1% = 1, 1%. = 1 and 
l-i = 2. Since for the other input parameters we find 
similar qualitative results, in order to avoid reputation, 
our analysis is restricted to the mentioned input param- 
eters. We find that variations of the mass of fragments 
with the ratio M Wjt otai/-M acc (i?disc) is more sensitive to 
the higher values of the exponent of viscosity, /i, and 
so we adopt // = 2 in our plots. Also, the ratio of gas 
pressure and total pressure, /3, is not very sensitive to 
the ratio of the total mass-loss rate and the accretion 
rate for different input parameters. The other input pa- 
rameters related to the outflow are changed to illustrate 
their possible effects on the physical properties of the 
system. The mass-loss power-law index £ is adopted 
values of 0, —0.5, —1 and —1.5 in Figures [T] and [5] and 
also, we have £ = — 2 and —2.5 in Figure [3] 

We also consider / = 1 and 5. As noted, radiation 
driven disc winds are expected to belong to the I = 1 
family of models and centrifugally driven disc winds are 
expected to belong to the I = 5 family of models. In 
our analysis, we will use the nondimensional factor Ze/e 
as a free parameter so that by changing this parame- 
ter we can consider appropriate values of the accretion 
rate. Evidently, higher the ratio Ze/e, larger the accre- 
tion rate. However, some authors introduce different 
values for the accretion rate. For example, for a central 
mass with mass M = 10 s M Q , Goodman & Tan (2004) 
proposed Ze/e = 10. Thus, in our analysis, the chosen 
values of Ze/e = 1 and 10 are acceptable, also we are 
using values of 0.03 and 0.3 for ao- 

Using equation ([30|l . we can calculate the mass of 
the first fragments for the set of the above input pa- 
rameters. Figure [1] shows the mass of the fragments at 
the self-gravitating radius versus the ratio of the total 
mass loss rate and the accretion rate for different in- 
put parameters. Each curve is represented by a pair 
of Z and the mass-loss power-low index £ as (Z,£). For 
a given set of the input parameters, a higher accretion 
rate implies fragments with higher mass. For exam- 
ple, while for a case with ao — 0.03 and Ze/e = 1, the 



mass of the first clumps will be between 6 x 10 2 M Q and 
1.4 x 10 3 M Q depending on mass-loss rate by outflow, we 
see that for ao = 0.03 and Ze/e = 10 the mass of frag- 
ments will increase up to a value between 4 x 1O 4 M0 
and 9 x 10 4 M Q . Also, for a fixed accretion rate, the 
mass of the first clumps decreases with increasing ao- 
In all plots of Figure [T] for Z = 5, we see that mass 
of fragments is more sensitive to the variations of the 
mass-loss power-law index £ in comparison to the so- 
lutions with 1 = 1. In other words, when outflows are 
centrifugally driven, the mass of the fragments highly 
depends on the mass-loss index £. But for radiation 
driven outflows this dependence is weak. Also, in the 
case of centrifugally driven outflows, for a fixed ratio 
-Mw,totai/-A^acc(-Rdisc) the mass of fragments increases 
with increasing the mass-loss index £. But in the case of 
radiation driven outflows we see an opposite variation, 
i.e. the mass of fragments decreases with increasing the 
mass-loss index £ if the other input parameters arc kept 
fixed. Also our plots show that for all input parame- 
ter when the ratio M Wit otai/-M acc (i?disc) increases, then 
mass of the first clumps decreases and interestingly the 
reduction is linearly proportional to the mass-loss rate, 
i.e. M frag — --4(Mw,totai/Macc(-Rdisc)) + B where con- 
stants A and B depend on the input parameters. 

Figure [2] shows the self-gravitating radius R sg (in 
Schwarzschild radius i?s) versus the ratio of the total 
mass loss rate and the accretion rate for different input 
parameters. For a given set of the input parameters, the 
self-gravitating radius increases with increasing the vis- 
cosity coefficient ao- The self- gravitating radius is not 
very sensitive to the variations of the mass-loss index £ 
when there are radiation driven outflows. When there 
are centrifugally driven outflows, the self-gravitating ra- 
dius increases with decreasing the mass-loss index, ex- 
cept for a case with small viscosity coefficient and high 
accretion rate (i.e., ao = 0.03 and Ze/e = 10). Also, the 
self-gravitating radius decreases with increasing mass- 
loss rate by the outflow unless the viscosity coefficient 
is small and the accretion rate is high. In the case of 
£ = —2 and —2.5, we can also determine self-gravitating 
radius and the mass of the fragments. Figure [3] shows 
the mass of the fragments (top) and the self-gravitating 
radius R sg (bottom) versus the ratio of the total mass 
loss rate and the accretion rate. Variations of the mass 
of the fragments and self-gravitating radius with the 
input parameters arc similar to Figures [T] and [2J 

4 Conclusion 

We presented a set of analytical solutions for the steady 
state structure of discs with outflows around supermas- 
sive black holes. Gravitational stability of the disc has 
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Fig. 3 The self-gravitating radius, R ss , (in Schwarzschild radius Rs) and the mass of the first clumps, Mf rag , (in solar 
mass) versus the ratio of the total mass-loss by the outflow to the accretion rate, M Wito tai/Afacc(-Rdisc), f° r M ~ 2.0, jix m = 0.6, 
k — 1, Ms = 1, / = \/2, rfo = f and ryk = 1. Here, the index of mass- loss rate by outflows is adopted as £ = —2 and —2.5. 
Solid and dashed curves are corresponding to the solutions with I = 1 and I = 5, respectively. Each curve is marked by a 
pair (Z,£). 
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been studied using our analytical solutions. We deter- 
mined the self-gravitating radius and the mass of the 
first clumps at this radius. We showed that as more 
mass, angular momentum and energy are extracted 
from the disc by the outflows, the mass of the first 
fragments decreases, though the self-gravitating radius 
increases except for a case with small viscosity coeffi- 
cient and high accretion rate. However, we think low 
values of cto are not acceptable in self-gravitating discs. 
We can conclude that the existence of outflows imply 
a more gravitationally stable accretion disc. Interest- 
ingly, there is a linear correlation between the mass 
of the fragments and mass-loss rate by the outflows, 
according to our solutions. The mentioned effects of 
the outflows on the gravitational stability of the discs 
should be considered in theoretical studies of star for- 
mation near to the supermassive black holes. 
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